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ABSTRACT 

We spectroscopically identified 7 massive evolved galaxies with magnitudes 17.8<K< 
18.4 at 1.3 < z < 1.7 over an area of ~ 160 arcmin 2 of the MUNICS survey. Their 
rest-frame K-band absolute magnitudes are -26.8<M^ < —26.1 (5L* <Lk <10L*) 
and the resulting stellar masses are in the range 3 -j- 6.5 x 10 11 M Q . The analysis 
we performed unambiguously shows the early-type nature of their spectra. The 7 
massive evolved galaxies account for a comoving density of (5.5 ± 2) x 1CP 5 Mpc~ 3 at 
(z) ~ 1.5, a factor 1.5 lower than the density ((8.4 ± 1) x 10~ 5 Mpc -3 ) of early-types 
with comparable masses at z = 0. The incompleteness (<~ 30%) of our spectroscopic 
observations accounts for this discrepancy. Thus, our data do not support a decrease 
of the comoving density of early-type galaxies with masses comparable to the most 
massive ones in the local Universe up to z ~ 1.7. This suggests that massive evolved 
galaxies do not play an important role in the evolution of the mass density outlined by 
recent surveys in this redshift range, evolution which instead has to be ascribed to the 
accretion of the stellar mass in late-type galaxies. Finally, the presence of such massive 
evolved galaxies at these redshifts suggests that the assembly of massive spheroids has 
taken place at z > 2 supporting a high efficiency in the accretion of the stellar mass 
in massive halos in the early Universe. 

Key words: Galaxies: evolution - Galaxies: elliptical and lenticular, cD - Galaxies: 
formation. 



1 INTRODUCTION 

The epoch of formation of high- mass (Mstar > W 11 Mq) 
early-type galaxies is one of the open questions relevant to 
the whole picture of galaxy formation and evolution. The 
uniform properties shown by the local ellipticals suggested 
the simple monolithic collapse scenario of galaxy formation 
(Eggen et al. 1962; Arimoto & Yoshii 1987). On the other 
hand, the recent picture outlined by hierarchical models 
(White & Frank 1991; Kauffmann 1996; Somerville & Pri- 
mack 1999) depicts the formation of local ellipticals through 
subsequent mergers: the higher the final stellar mass of the 
galaxy, the later it has been assembled. The most massive of 
them (10 n -10 12 M Q ) populating the brightest end (L>L*) 
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of the luminosity function of galaxies are those reaching their 
final mass most recently in this scenario, possibly at z < 1. 
Therefore, looking for z > 1 early-types with stellar masses 
comparable to those of the most massive local ones is one of 
the most direct ways to address the question of galaxy for- 
mation. This is what we are doing through a near-IR spec- 
troscopic survey of early-type galaxy candidates selected to 
be at z > 1 and to have stellar masses well in excess to 10 
Mq (Saracco et al. 2003a,b). The candidates consist of a 
complete sample of 31 bright (K'< 18.5) Extremely Red Ob- 
jects (EROs) with colours R-K'^ 5.3 selected over two fields 
(~320 arcmin 2 ) of the Munich Near-IR Cluster Survey (MU- 
NICS; Drory et al. 2001). Here, we report the spectroscopic 
confirmation for a sample of 7 high-mass {Mstar > 3 x 10 
Mq) field early-type galaxies identified at 1.3 < z < 1.7. 
Throughout this paper we assume Ho=70 km s _1 Mpc -1 , 
Qo = 0.3 and A = 0.7. 
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Figure 1. Expected K-band apparent magnitude for an 
\jk =1.5L* (thin line) elliptical and for a Ljf =2.5xL* (thick line) 
elliptical as a function of rcdshift. The elliptical has been modeled 
by assuming a SSP model based on the evolutionary population 
synthesis of Maraston (1998, 2004). The elliptical is 2 Gyr old at 
z ~ 1.5 and evolves passively down to z = 0. The dashed line 
marks the selection criterion K^18.5. 



2 SAMPLE SELECTION, OBSERVATIONS 
AND SPECTROSCOPIC CLASSIFICATION 

2.1 Sample selection 

The main aim of our spectroscopic survey is to identify 
early- type galaxies at z > 1.2 having stellar masses com- 
parable to the most massive early-types in the local Uni- 
verse (10 11 — 10 12 Mq). The red optical-to-near-IR colour 
(R-K^ 5.3) favors the selection of z > 1 passively evolved 
galaxies while the magnitude K'< 18.5 assures the selection 
of systems with stellar masses well in excess of M sta r > 10 
Mq. Indeed, pushing back in time a local massive elliptical 
(L=2.5L*, 2 — 3x 10 11 M Q ) it would be observed with a mag- 
nitude brighter than K=18.5 up to z ~ 1.5. This is shown 
in Fig. 1 where the expected K-band apparent magnitude of 
local bright ellipticals is plotted as a function of redshift (we 
considered M^=-24.3, Kochanek et al. 2001) in case of pure 
passive evolution. We modeled the elliptical by assuming 
a Simple Stellar Population (SSP) model (Maraston 1998, 
2004) 2 Gyr old (Z=Z0) at z ~ 1.5. It is worth of noting that 
even little star formation would produce a brightening of the 
expected apparent magnitude. Thus, K'<18.5 is a stringent 
criterion which, together with the colour criterion, assures 
that the resulting early- types at 1.2 < z < 2 are as massive 
as the most massive ones (10 11 — 10 12 Mq) in the present- 
day Universe. The resulting sample of EROs comprises 31 
galaxies. Thus, we measured a surface density of R-K^5.3 
galaxies of 0.10±0.02 arcmin -2 at K< 18.5, in agreement 
with previous estimates over larger areas (e.g. Daddi et al. 
2000; Martini 2001). 

2.2 Observations 

The spectroscopic observations, with typical exposure of 
about 4 hours for each source, were carried out in Octo- 



ber 2002 and November 2003 at the Italian 3.6 m Tele- 
scopio Nazionale Galileo (TNG). The prism disperser Am- 
ici mounted at the near-IR camera NICS of the TNG was 
adopted to carry out the observations. This prism provides 
the spectrum from 0.85 fj,m to 2.4 /jm in a single shot at a 
nearly constant resolution of A/AA ~ 35 (1.5" slit width). 
This resolution is best suited to describe the spectral shape 
of sources over a wavelength range of ~ 15000 A and to de- 
tect strong continuum features such as the 4000 A break 
in old stellar systems at z > 1.2. On the other hand, the 
extremely low resolution makes unfeasible the detection of 
emission/absorption lines for sources as faint as our galax- 
ies. Until now, we carried out spectroscopic observations for 
~ 60% of the whole sample identifying 10 early-types. The 
analysis of the spectral properties of the whole sample of 
early-types is presented in a forthcoming paper (Longhetti 
et al. 2004). Seven of the 10 early- types fall on one of the two 
selected fields, the S2F1 field (~ 160 arcmin 2 ), where we col- 
lected 13 spectra out of the 19 EROs satisfying the selection 
criteria. In Tab. 1 we report the broad-band photometry of 
the 7 early- types in this field. 

2.3 Spectral Classification 

In Fig. 2 the (smoothed) near-IR spectra (black thick his- 
togram) of the 7 galaxies are shown. The spectra drops very 
rapidly at ~ 0.9 — 1 fim concurrent with the 4000A break 
placing the galaxies at z > 1.3 and suggesting an early- type 
spectral nature. We searched for the best-fitting template by 
comparing the observed SED of each galaxy, constituted by 
the broad band photometry (B, V, R, I, J and K) and by the 
observed near-IR continuum (from 0.9 /im to 2.3 /im), with 
a set of spectrophotometry templates (Z=Z0 and Salpeter 
IMF) based on the Bruzual & Chariot (2003) models and 
on the library of Simple Stellar Population (SSP) of Maras- 
ton (1998, 2004). The best-fitting templates have been ob- 
tained through a \ 2 minimization applied over the whole 
wavelength range. Besides the SSPs, models with declining 
star formation with time scales r in the range 0.1-4 Gyr 
has been considered. All the observed SEDs are best-fitted 
by SFHs with very short time scales (r ^ 0.3) Gyr and 
extinctions E(B-V)< 0.15 (with the Calzetti et al. (2000) 
extinction law), providing ages in the range 1.5-4 Gyr. A 
comprehensive study of the properties of the stellar popu- 
lations in our massive early-types and of their dependence 
on different IMFs and metallicity is presented in Longhetti 
et al. (2004). In Table 1, the spectroscopic redshift and the 
best-fitting SFHs are summarized. It is worth of noting that 
a set of SSPs with ages in the range 1-3.5 Gyr provides a 
good fit as well to all the observed SEDs. To rule out the 
possibility that some (or all) of the observed SEDs can be 
fitted by a young dusty starburst, we forced the fitting pro- 
cedure to search for an acceptable fit among a set of tem- 
plates made up by the 6 starburst templates (SB1-SB6) of 
Kinney et al. (1996) and a starburst model described by 
a constant star formation rate (est). Extinction has been 
allowed to vary in the range 0<E(B-V)< 2 in the fitting 
procedure. We did not obtain an acceptable fit for any of 
the 7 galaxies. In particular, we fail in simultaneously fit- 
ting the red part of the spectrum at A > 1.4 fim and the 
blue part at A < 0.9 fim. This is shown in Fig. 2, where we 
plot the best-fitting starburst template (dotted line) to the 
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Figure 2. Near-IR spectra (black histogram) of 4 out of the 7 early-type galaxies. The mean observed spectrum of local ellipticals 
(thin grey line) of Mannucci et al. (2001) and of Coleman et al. (1980) are superimposed on the observed spectrum of S2F1_357 and of 
S2F1_633 and on the observed spectrum of S2F1_443 respectively. A SSP 1 Gyr old at solar metallicity is superimposed on the observed 
spectrum of S2F1_511. For the other galaxies, we superimposed on the Amici spectra the best-fitting template (see Tab.l). In the case 
of S2F1_527, the best-fitting starburst template is also shown (dotted line). The horizontal errorbars represent the atmospheric windows 
with an opacity larger than 80%. The filled symbols are the photometric data in the B, V, R, I, J and K' bands from the MUNICS 
catalog (Drory et al. 2001). 



most favorable case represented by S2F1_527. Thus, dusty 
starbursts do not reproduce the sharp and deep drop seen 
in the continuum at 0.85< A i, s < 1.0 /xm. This confirms the 
early-type spectral nature of the 7 galaxies. For comparison, 
in Fig. 2, the mean observed spectrum of local ellipticals of 
Mannucci et al. (2001) (thin grey line) is superimposed on 
the Amici spectrum of S2F1_357 and of S2F1_663, while that 
of Coleman et al. (1980) is superimposed on the spectrum of 
S2F1_443. For the other galaxies, the best-fitting template 
is superimposed on the observed spectrum. It is worth of 
noting that, the highest redshift evolved galaxy of our sam- 
ple (S2F1_443) is also one of the X-ray emitting EROs we 
detect on the S2F1 field (Severgnini et al. 2004). 



IMF) and by the scale factor applied to the flux of the red- 
shifted template to fit the observed fluxes. This latter esti- 
mate does not consider a certain band but takes into account 
the whole observed SED. The two estimates provides sim- 
ilar values (within few percent) since most of our data are 
in the near-IR and the uncertainty on them are lower than 
the uncertainty affecting the optical data. Consequently, the 
fitting procedure weights mostly the near-IR data. An anal- 
ysis of the dependence of the stellar mass of the 7 evolved 
galaxies on the IMF and on the metallicity is discussed in 
Longhetti et al. (2004). In Tab.l we report for each galaxy 
the rest-frame K-band luminosity, the M/Lk and the rele- 
vant stellar mass. As expected, all the galaxies have stellar 
masses in the range 3 -v- 6.5 x 10 11 M©. 



3 LUMINOSITIES AND STELLAR MASSES 

The bright K'-band magnitudes (17.8<K'< 18.4) of our el- 
lipticals and their redshift imply luminosities L^>L*. Since 
local L* galaxies have stellar masses of the order of 10 11 
Mq, we expect all the 7 early-types more massive than 10 11 
Mq leaving aside any model assumption. In order to derive 
the rest-frame K-band luminosity, the proper k-correction 
of each galaxy has been computed by means of the rele- 
vant best-fitting template. Each best-fitting template has 
been multiplied with the transmission curve of the K filter 
to derive the k-corrections. In Tab. 1 we report the K-band 
absolute magnitudes thus obtained. As expected, the 7 early- 
types have rest-frame near-IR luminosities 5L* <Lk ^10L*. 
The stellar mass M s tar of each galaxy has been derived by 
the K-band luminosity through the mass-to-K-band light ra- 
tio relevant to the best-fitting template (Z=Zq and Salpeter 



4 THE DENSITY OF VERY MASSIVE 
EARLY-TYPE GALAXIES 

We estimated the co-moving spatial density of the 7 massive 
early-type galaxies and its statistical uncertainty as: 

o=£^ ^HE(t^) 2 F /2 (1) 

* max L r max -I 

i i 

where 

Vmax = -J— I ^-d,Z (2) 

47T / dz 

J zi 

is the comoving volume. The solid angle ui subtended by 
the S2F1 field is ~ 1.3 ■ 10 -5 strd and z m ax is the maxi- 
mum redshift at which each galaxy would be still included 
in the sample. In the derivation of the z m ax of each galaxy 
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Table 1. Properties of the 7 early-type galaxies. Magnitudes are in the Vega system. The SFHs, the stellar masses and the mass-to-K- 
band light ratio refer to Z=Zq models with Salpeter IMF. 



ID 


K 


R-K 


Zspec 


M K 


SFH 


■Mstar 


M/L K 


(L*=o/L*) 












[Gyr] 


[10 11 M ] 






S2F1.357 


17.84±0.08 


6.0 


1.34±0.05 


-26.6±0.12 


SSP 


5.0 


0.5 


3.0 


S2F1.527 


18.30±0.15 


> 5.7 


1.35±0.05 


-26.3±0.20 


T = 0.1 


3.0 


0.4 


2.3 


S2F1.389 


18.23±0.12 


5.5 


1.40±0.05 


-26.5±0.15 


t = 0.3 


3.5 


0.4 


2.7 


S2F1.511 


18.14±0.15 


6.1 


1.40±0.05 


-26.2±0.20 


T = 0.1 


3.0 


0.4 


2.1 


S2F1.142 


17.84±0.07 


6.0 


1.43±0.05 


-26.6±0.12 


r = 0.3 


6.5 


0.6 


3.0 


S2F1.633 


18.20±0.12 


> 5.7 


1.45±0.05 


-26.1±0.15 


T = 0.1 


4.0 


0.6 


1.9 


S2F1.443 


18.40±0.15 


> 5.6 


1.70±0.05 


-26.8±0.20 


T = 0.1 


5.0 


0.4 


3.6 



we computed the k-correction by using the relevant best fit- 
ting template. The lower bound in the integration is set to 
21 = 1.2. It is imposed by the adopted spectroscopic wave- 
length range which does not allow us to detect the Balmer 
break of galaxies at z < 1.2. We find that, at the aver- 
age redshift (z max ) ~ 1.55, the 7 massive early-types ac- 
count for a co-moving density p — (5.5 ± 2) x 1CP 5 Mpc -3 
over a volume of about 1.5 x 10 5 Mpc 3 . Given the stellar 
masses of the 7 early-types the resulting mass density is 
pstar = (2.3 ± 0.9) x 10 7 M Q Mpc" 3 . Such densities are 
likely to be lower limits because of the incompleteness of 
our spectroscopic observations (13 spectra collected out of 
the 19 EROs in the field). In order to compare these den- 
sities with the local values we integrated the K-band LF of 
local early-type galaxies derived by Kochanek et al. (2001) 
described by M^ = -24.3 ± 0.06 and $* = 1.5 ± 0.2 x 10~ 3 
Mpc -3 . The lower bound to the integration has been de- 
rived in two independent ways. The first one lies in deriving 
the lowest luminosity that our galaxies would have at z=0. 
This has been obtained by assuming a pure passive evolu- 
tion from z apec to z = 0. The luminosities thus obtained 
are reported in Tab. 1 and place the lower limit at 1.9L* as 
imposed by S2F1_633. In the second way, we derived the K- 
band luminosity from the lowest stellar mass of our galaxies 
(3 x 10 11 Mg) by means of the maximum mass-to-light ra- 
tio that our galaxies could have at z — 0. Since their stellar 
populations have to be formed at z > 2, they would be ~ 11 
Gyr old at z = 0. Thus, the maximum mass-to-light ratio is 
M/L K ~1.6, given by a SSP 11 Gyr old with Salpeter IMF. 
Also in this case we obtained a luminosity of about 1.9L*. 
It is worth of noting that other IMFs would produce lower 
values of M/Lk (e.g. M/Lk — I-S with Kroupa IMF and 
M/Lk=1-1 with Miller-Scalo) and, correspondingly, higher 
luminosities (2.3L* and 2.7L* respectively). By integrating 
the local LF of early-type galaxies at luminosities brighter 
than 1.9L* we obtained a density of (8.4± 1) x 10~ 5 Mpc -3 . 
We summarized our results in Fig. 3 where we report also 
the comoving density we derived for the 3 early-types with 
L z= o > 3L*. The density we estimated at z ~ 1.5 repre- 
sents 65% of the density of their counterparts at z = 0. Even 
if this difference is not statistically significant, this suggests 
that three more massive early-types should be expected over 
the 160 arcmin 2 in case of no density evolution. On the other 
hand, our estimate is affected by a spectroscopic redshift in- 
completeness of about 30%. When this incompleteness is 
taken into account the observed difference tends to vanish. 
Thus, we conclude that the number density of early-types 



with stellar masses comparable to the most massive early- 
types populating the local Universe does not show evidences 
of decrease up to z ~ 1.7. 

In Fig. 3 we also report the density relevant to the four 
old spheroidal galaxies (empty circle) spectroscopically iden- 
tified at z > 1.6 by Cimatti et al. (2004), the density of mas- 
sive (10 11 M©) evolved galaxy candidates at z > 2 (open tri- 
angle) found by Saracco et al. (2004) on the HDF-S and the 
density found by Caputi et al. (2004) on the GOODS/CDF-S 
area (crosses), both on the basis of photometric analysis. Fi- 
nally, the upper limit to the number density of galaxies with 
masses larger than 2 x 10 11 Mq (open squares) by Drory et 
al. (2004) is shown. 



5 CONCLUSIONS 

We spectroscopically identified 7 bright (17.8 <K< 18.4) 
massive evolved galaxies at < z < 1.7 over an area of about 
160 arcmin 2 of the MUNICS survey. These galaxies turned 
out to have rest-frame K-band luminosities 5.5L<Lk ^11L* 
and stellar masses in the range 3 — 6.5 x 10 11 Mq. At the 
mean redshift of z ~ 1.5 these 7 early-types sample a vol- 
ume of about 1.5 x 10 5 Mpc 3 and account for a comoving 
number density p = (5.5 ± 2) x 10~ 5 Mpc~ 3 and a stellar 
mass density p sta r = (2.3 ± 0.9) x 10 7 M Mpc" 3 . These 
densities represent 65% of the values at z — for early-types 
with comparable mass. The incompleteness of our spectro- 
scopic observations (30%) accounts for this deficiency. Thus, 
our results show that the number density of the most mas- 
sive early-types in the present-day Universe keeps essentially 
constant down to z ~ 1.7. This suggests that, massive early- 
types do not take part in the evolution of the stellar mass 
density in the redshift range < z < 1.7 and that the de- 
crease of the stellar mass density detected in this redshift 
range by recent surveys (e.g. Rudnick et al. 2003; Drory 
ct al. 2004; Fontana et al. 2004) has to be ascribed to the 
accretion of the stellar mass in massive late-type galaxies. 
This qualitatively agrees with the concurrent increase of the 
cosmic star formation rate due to late-type galaxies seen 
in the same redshift range (e.g. Madau et al. 1998). The 
high stellar masses and the number density of the 7 massive 
evolved galaxies imply that they were fully assembled at the 
observed redshift pushing their formation at z > 2, as also 
suggested by the recent discovery of an evolved spheroidal 
galaxy at z ~ 1.9 (Cimatti et al. 2004). This is in agreement 
with the results based on the spectral analysis of the stellar 
populations in our sample of massive early-types at z ~ 1.5 
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Figure 3. Number density of galaxies as a function of redshift. 
The large filled symbols represent the densities we derived at z ~ 
1.5 for the 7 early-types with L z — o > 2L* (circle) and for the 
3 with L z= o > 3L* (square). The small filled symbols at z = 
represent the number density of E/S0 galaxies brighter than 2L* 
(circle) and 3L* (square) respectively obtained by integrating the 
local LF of galaxies of Kochanek et al. (2001). The open circle 
is the density of spheroidal galaxies by Cimatti et al. (2004), the 
open triangle and the crosses are the density of M. > 10 11 Mg 
evolved galaxy candidates by Saracco et al. (2004) and by Caputi 
et al. (2004) respectively and, finally, the open squares are the 
upper limit to the number density of galaxies with M > 2 X 10 
M by Drory et al. (2004). 



which suggest a formation redshift z/ > 2 (Longhetti et 
al. 2004). It is worth of noting that similar results are also 
derived by the analysis of the absorption lines indices of lo- 
cal samples (Thomas et al. 2002, 2004) which suggest short 
timescales (~ 0.4 Gyr) and high-z of formation for mas- 
sive early-types. Thus, at variance with the expectations of 
hierarchical models, the most massive early-types in the lo- 
cal Universe do not seem to be the last galaxies to complete 
their assembly. The high formation redshift suggests an high 
efficiency in the accretion of the stellar mass of early-types in 
the early Universe and an high star formation preferentially 
in massive haloes. 
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